TITLE: 

TYPE: 

NUMBER OF WORDS : 

TEMPORARY STORAGE: 

ACCURACY: 

DURATION: 

DESCRIPTION: 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

1/2 sin X/X (SADOI Only) 
Closed subroutine (Uses T 5) 

37 

0, 1, 2, and 3. 
-37 



LIBRARY ROUTINE T 7 - 293 



+ 2 

(15 + |) milliseconds for |x| > l/2, where 2~ n |x| = 1. 

7 milliseconds for |x| < l/2. 

X is treated as a double precision number with an integer 

and fractional part. Enter with the integer part of X 

in location 3 and the fractional part of X in the accumulator, 

The accumulator will contain 

1/2 53fJL 

upon exit. 

This routine uses library routine T 5 - 157 and identifies 
it by means of the symbolic address (T5). (T5) must be 
supplied by the programmer 

For |x| > 1/2, T 5 is used to compute 1/2 sin X. This 
quantity is shifted right simultaneously with X until 
division by X is possible. For X < 1/2, a power series 
is used. The power series coefficients were obtained 
using routine KA1. 
NUMBER REPRESENTATION: 

The double precision numbers are correctly represented by 
a signed integer and a positive fraction. For example, 
the number +2.4 should be represented by the integer 
+2 and the fraction +0.4. The number -2.k should be 
represented by the integer -3 and the fraction +0.6. 
Numbers generated by multiplication are correctly re- 
presented since q^ = after multiplication in all cases. 



CAUTION: 



METHOD: 



DAT E March 30, i960 
PROGRAMMED B Y^.Z^MMjl 



APPROVED BY 



ns 



f 3 ^ 






LOCATION 


■ ■ ■ ■ - — 
ORDER 


— 


i 
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; 




00K(T7) 






: 





4o 36L 
K5 31L 




save fraction 




1 


42 20L 
L5 L 




link 




2 


^2 24L 










L5 3F 




integer part of X 




3 


32 6L 
LI 36L 




) 




4 


10 39F 
LO 3F 




I convert negative 
[ X to positive 




5 


40 3F 
01 39F 




/ 




6 


4o 36L 








43 18L 








7 


L3 3F 
36 21L 




integer part = 0? 




8 


5L 36L 
7J 30L 




order pair near 3/4 




9 


50 3F 








10 


74 30L 
00 39F 
50 10L 




i[x l + 2 -39 XF ] 
into A. 




11 


26 (T5) 

40 IF 




1/2 sin X in 1. 




12 


50 36L 
F5 18L 


■ 


*} 




13 


42 18L 
L3 3F 




/ Shift X until 




14 


32 i6l 
L5 3F 




' pure fraction. 
Count shifts. 




15 


10 IF 
40 3F 


V 


/ 




16 


22 12L 








S5 F 
— 


from 14L 







LOCATION 


ORDER 




NOTES PAGE 2 


T 7 


17 


40 3F 
50 19L 




save shifted X. 
Clear Q. 




18 


L5 IF 
10 F 


» 


I Shift 1/2 sin X same as X. 




19 


00 IF 
66 5F 


1 


) 




20 


S5 F 
22 F 




exit . 




21 


50 56L 
79 56L 


from 7L 






22 


40 F 
LO 8l 




-X/ in 




25 


56 8L 
50 5F 




jump if X > 1/2. 
J 5 = initially 




24 


79 F 
L4 F 




1 




25 


40 5F 

TT>C o)iT 

.r ^ t-ru 




Compute l/2 y 
r^„ 




26 


42 24L 
LO 29L 




series when 




27 


52 25L 
L5 5F 




|X < 1/2 




28 


22 20L 
00 F 








29 


L9 F 
L4 56L 




end test 
constant 




50 


00 F 










00 318 509 886 184 J 


1 




51 


00 F | 

00 1 570 058 J 






52 


80 F | 








00 999 900 795 557 J 


(coefficients 




55 


00 F J 


[for power 






00 4 166 666 514 j 


li 





LOCATION 
36 



ORDER 

80 F 

00 916 666 666 671 J 

ho F 

00 F 

00 F 

00 F 



NOTES 
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T 7 



series. 



o- 



Fractional part of X, 



